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Abstract 

The aim of this paper is to study some examples of exponentially 
harmonic maps. We study such maps firstly on flat euclidean and 
Minkowski spaces and secondly on Friedmann-Lemaitre universes. We 
also consider some new models of exponentially harmonic maps which 
are coupled with gravity which happen to be based on a generalization 
of the lagrangian for bosonic strings coupled with dilatonic field. 



1 Introduction 

Exponentially harmonic maps were introduced by James Eells and studied by 
J. Eells and L. Lemaire 0. These maps generalize usual harmonic maps || 
in the following sense. Let (M, g) and (N, h) be two Riemannian manifolds 
and <f) : M — > N : x — > 4>{x), a smooth map. An exponentially harmonic map 
is then an extremal of the following functional: 

E(<f>) = [ exp(e(0))d/i(0) (1) 

J M 
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where dfi(<j)) is the riemannian volume element and 



Lb 



is the so called "energy density" of the map <fi. In local coordinates and 
(4> a ), the equation for exponentially harmonic maps, which derives from the 
variation of the functional flU), is 



( U\ J ^ / d 2( ^ (M) 90 a (iV) «90 fe d<f) c \ 

+a^h h T b ^^^-^-^-\-0 (3) 
+9 9 flfcl* q^q^q^q^]-^ ^ 

where r^Jg and T^" 1 are the Christophel symbols of the Levi-Civita connec- 
tion on M and N. This equation involves, as a particular case, the equation 
of usual harmonic maps. If we drop the exp(e(0)) factor and restrict ourselves 
to the first three terms we get: 

y \dx a dxP a(3 dx-r bc dx a dx?) K > 

which is the field equation for a usual harmonic map (which is nothing but 
a non-linear sigma model in the physicist's language). When N = R, (^) is 
simply written: 



1 d 



V99 a ^ ) = (5) 



which is the Laplace equation in local coordinates. Following the work of 
Luc Lemaire, it is important to note, however, that the properties of the 
exponentially harmonic maps are very different of those of usual harmonic 
maps. This comes from the fact that the functional ([!]) changes completely 
when we perform a conformal shift on the metric h: h — > X.h. Then ([!]) 
becomes: 



E\{4>) = / [exp(e(0))] d[i((f)) X = constant. (6) 

Jm 

In the case of usual harmonic maps this metric shift has no influence on the 
harmonic maps equations, the functional is simply multiplied by a constant. 
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In the sequel we will study somme particular cases of (|3]) even when M 
is noncompact and (jH) unbounded. In this case, equation (^) is taken as the 
definition of what we call an exponentially harmonic map. Following Eells 
and Lemaire, we will consider the energy-momentum tensor asociated to (f>: 

T M) = ex p( e (^)) (9^ - g^a^A fe) ( 7 ) 

which is conserved (5 ,a/3 V Q T a ^(0) = 0) when is a solution of (|3|). 

It is worth noting that some particular cases of (|3|) were studied by mathe- 
maticians in the context of the theory of elliptic partial differential equations. 
Gilbarg and Trudinger [[| quote, for example, the equation: 



dx l dxi dx l dxi 
m'=i 



deriving from the variational principle: 8 J" exp(| ^^(J^) 2 )^ 1 ■■■dx n 

Lut us remark finally that, if we want to make some contact with physics, 
we have to modify (|I|) as follows: 



E' x = / (exp(Ae(0)) - 1)^(0) 

J M 

1 ■ A 1 • • .. > 



A 



dn(<l>) ( 9 ) 



' M L 

which, if we interpret (<f) a ) as a spin vector, looks like the Hamiltonian for 
the continuum limit of a Heisenberg model. When A is small enough, the 
variational principle 5E' X = leads to equations which approximate those of 
the usual harmonic maps (sigma models). The equations derived from E' x 
can be obtained directly from via the shift cj) -»• y/\(f) (A > 0). 



2 Exponential harmonic maps on flat spaces 

In order to get some feeling about the solutions of (f|), we study here scalar ex- 
ponentially harmonic maps on two-dimensional Euclidean and Minkowskian 
manifolds: 

: E -»• R E = R 2 or R 1 ' 1 . 

2.1 E = R 2 

Let = (x,y), <f) xf i = dfp/dx^ and 4> x n x v = d 2 (f>/dx^dx u . Then, equation 
(|3|) becomes 

(1 + <f>l)<f>XX + 2</> x <f> v <f> SBy + (1 + <pl)<Pyy = 0. (10) 
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We are searching solutions of the following form: (f>(x, y) = F(x) + G(y). We 
are thus led to the equations: 

[1 + {F x f\F xx = -[1 + (G y ) 2 ]G yy = I (11) 

where I happens to be an arbitrary positive real constant. I = leads, of 
course, to trivial solutions y) = ax + by. If / > 0, let p = F x . This leads 
to the system: 



F x = p(x) (12) 
I 

1 + p 2 



p x = -r— r (13) 



Then (O) gives 



x = -j{-p 3 + p — c); c = constant. (14) 
Let us derive ( |i~4"D with respect to F. This gives (with fc=constant): 

% = Y x = l {l+T>2) % and P 4 + V-4/(F + A;) = 0. (15) 
The elimination of p from flT4]) and (|15D gives: 

F(x) = ^{[S + (x;l;c) + S4x;l;c)] 4 

+2[S + (x;l;c) + S_(x; Z; c)] 2 } -k (16) 
p = S + (x;l;c) + S-(x;l;c). (17) 

where S±(x; I; c) = {|(c + Ix) ± (1 + |(c + Za;) 2 ) 1 / 2 } 1 ^ 3 . Similary we get the 
solution G(y). A solution of ( |T0D can thus be written, with the above nota- 
tions as, 

G(y) = -^{[S + (x;-l;c) + S4x;-l;c)]' 

+2[S + (x; -I; c) + S.(x; -I; c)] 2 } - k 

which leads to the solution we look for. 

If we write q(y) = S + (x; —I; c)+S-(x; —I; c), then we can represent 0(x, y) 
in the following parametric form: 

x= y(y + P-c), V = — (y + ff-c), 
0(x, y) = ^(P 4 + 2p 2 -q 4 - 2q 2 ) + est. 



It is also possible to solve ( |I0D by the so-called hodograph method ||. Let 
u = (p x ,v = 4> y and x = x(u,v),y = y(u,v). If J = x u y v - x v y u ^ 0, we get 
immediatly: 

u x = y v /J, v x = -y u /J 

Ziy Xyj Vy X U j J 

We know that if v x = u y then x v = y u . Thus, there exists a function f(u,v) 
such that x = f u and y — f v . Equation (|i"0[ ) can be written 



(1 + u )f m - 2uvf {uv) + (l + v )f uv = 
or, using polar coordinates, u = r cos 6,v = r sin 8: 



frr + (r + -)fr + (l + \)fee = (19) 



which can be solved by factorizing: f(r,8) = R(r)T(8). We have now to 
solve the two equations (with an arbitrary constant a): 



R rr + (r + -)R r - a 2 (l + — )R = (20) 

T ee + a 2 T = 0. (21) 

If we are interested in periodic solutions, we set T{6) = A cos(a9) +B sin(a6 l ). 
Thus (p0|) can be reduced to the normal form by the following functional 
change: 

1 r 2 ~ 

R^R; R(r) = — exp( )R(r). 

i r 4 



Equation (|2Cf ) then becomes: 



1 - 4a 2 1 



r 



Rrr + + a 2 ) + (— |— )^ -j) R - (22) 

This equation can be reduced to the Whittaker equation 0. Writting 
R(r) = r -1 / 2 M(r 2 /2), then the function M(e) sastifies the Whittaker equa- 
tion: 

d 2 -1 a 1 - u 2 

M(e) = M a ^ /2 (e); -^M a ^, 2 + (_ + - + __£_)M aiAt/2 = 

with a = ±|(a 2 + 1), n = ±a. This solution can be expressed in terms of 
the confluent hypergeometric function xF\ which leads to the solution of (|20|) 
written as 

R (r) = r iFl{l + 2 + T ' 1 + a ' r /2) - 

5 



If a > 0, this solution is regular at the origin. We know that iiq(m, m, x) = 
exp(x). Then, for a = 1, we get R(r) = r/2, which is a solution of (f20l). 
Nevertheless this trivial solution (and related to the case a = 0) leads to 
some problems because, if we set T{6) = cos 9, f(r, 6) = r cos 9, then /„ = 0, 
J = 0. 



2.2 E 



.1,1 



The computations are similar to the preceding case. The equation @ leads 
to 

(1 + (f>l)(j) xx - 2(j) x (f)y(f) x y - (1 - <f>l)(j)yy = 0. 

The solutions of the form <f)(x,y) = F(x) + F(y) can be written with a 
parametric representation: 

!/P 3 , \ 1 /-9 3 , n 

^ = ^(y+P-c), 2/ = -r(-g- + g-c), 

0(x, y) = ^-(p 4 + 2p 2 - g 4 + 2g 2 ) + est. 

3 Exponentially harmonic maps on Friedmann- 
Lemaitre universe 

Let M be a Friedmann-Lemaitre (FL) universe endowed with the following 
metric: 

ds 2 = dt 2 - R 2 (t) ( dr ] - + r 2 (d6 2 + sin 2 8d(J) 2 )) . (23) 
\ 1 — kv J 

We consider exponentially harmonic maps: : M — > R which are, for physi- 
cal reasons explained above, extremals of the functional E^. They satisfy a 
modified version of equation flU) which is written in this case: 

0(1 + A0 2 ) + 3^ = 0, (0 = t , i? = i? t ) (24) 
if we restrict ourselves to — 4>{t)- This gives: 

R 3 (t) = -^-exp( 2 ), (a : positive constant). (25) 

|0| 2 
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Let us take, for example, R(t) = R {j-) 2 ^ i.e., M is Euclidean FL-universe 
(without cosmological constant) and t Q = 2/3h Q where H Q is the present value 
of the Hubble constant (R/R\ t=to = H Q ). We get the field equation 

|0|exp(^0 2 ) = l/&t 2 , (b = R 3 Jat 2 ). 

i) When t — > oo: this gives 0(i) — > constant. 

ii) When t —>■ and A small: 0(t) ~ O ± l/&£ which is not regular at t — 0. 
See Figure 1 

When A ~ we recover the well-know theory of harmonic maps coupled 
with gravity studied by T. Hirschmann, R. Schimming ||, T. Hughes, T. 
Kato, J. Marden and Lemaire and A. J. Vanderwinden [|J. 

4 Exponentially harmonic maps on F-L uni- 
verse coupled with gravity 

Let us start now with the following action, describing a gravitational field 
coupled with an exponentially scalar field: 



5(0) = ~J V=gd A x Mr- exp(^ a 09» - Aj + c mat 



(26) 



mat) 



K is a constant, A is the cosmological constant, C mat is the Lagrangian 
density for matter and <fi : M — > K. is a scalar field defined on a 4-dimensional 
spacetime. The variation of S(<j)) leads to Einstein's equations: 

" \R9»v = \ (-e^ a * - A) - A^0e^^} + T* 

where T^ at ^ is the usual energy-momentum tensor for matter. The variation 
with respect to <fi gives a field equation which is very similar to (|3|). Let us 
assume that <fi = <fi(t) and let M be a Friedmann-Lemaitre universe. Then 
the field equations can be written: 

3i*v+4 =A -M e ^ (1 _ A ^_f) <27> 

R k ( 1 A 5 2 A\ 

0(1 + A0 2 )+ 3^0 = 0. (29) 
R 

Let us consider some particular cases. 
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4.1 A flat F-L universe without matter 



0(1 + A0 2 )+ 3^0 = 0. 

Let us define y = 4> an d let H = be the Hubble constant. We check that 
H + H 2 = J|, then the equations above define a dynamical system: 



H = -^y\^ (30) 

y = ~TTv (31) 



submitted to the constraint: 



# 2 = f ((-1 + Ay 2 )e^ 2 - A) . (32) 

Let us derive ([32]) and, using (|3~lD , let us compare the result with (|^). We 
see that the preceding equations are compatible. We are thus led to solve 
the following equation: 



? (1 + Ay 2 ) - ( V -) (1 - \y 2 ) = -K\y\hv\ (33) 

y \yj 4 

For very small values of A, this becomes: 

i-(i) a *0. (34) 

y \yJ 

This gives solutions: 

(j)(t) = ce at + e; R(t) = be Ht 

where a, c are constants which the same sign, b a positive constant and e is 
an arbitrary constant. H — — | (H is a first integral when A = 0). The 
equation ( |3"2"| ) allows to write A = — 1 — See Figure 2 
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It is interesting to note here that the coupling of with the gravitational 
field can make </> regular at t = in this Euclidean case, which was not the 
case in the uncoupled situation (see section3). 

Let us return to the case where A is an arbitrary constant. We define 
z — y/y, equation (|33"D leads to 



y 

z 



2Xy 2 ( 3K a 2 2 



l + Ay 2V 8 

which is a dynamical system with a fixed point (y,z) = (0,0). If we put 
u = y 2 we get: 



(35) 

1 + ^8- ^ (36) 
The Hubble constant can be deduced from the numerical integration of this 



u 

z 



2zu 

2Xu ,3K a, 
( e 2 



dynamical system. From (R2) we get 



H = ± 



K 
~6 



((1 - Xu)e^ u + A) 



1/2 



(37) 



Around u ~ (i.e. = constant) we can neglect terms of order greater or 
equal to 2. Using ( |35l) - (|36|) we have at first order: 



««-^» (|A|(-^ + 2^ 2 ))+ Mo . 



(3? 



As u = y 2 , we have to suppose that (for A > 0): — ^f- + 2z 2 < j- when u Q = 0. 



4.2 A curved F-L universe without matter 

Let us derive (|29"D. Using ( p7|) and (p8|) with y = 0,we find the following 
equation: 

2 



12(1 + A 2/ 2 )^ + 4(A 2/ 2 -AV- 4)(2- 



V 

y ' \y, 

-3K(Xy 2 + 2)e^ 2 - 6KA = 0. 



(39) 



Equation (|39| ) can be solved numerically using the variable change: u = y 2 . 
Then (B9f) becomes: 



6(l + Au)--(AV-2Au + 7) - - 3K(Xu + 2)e* u - QKK = 0. (40) 
it V it 



Making the substitution z 

2zu 



lit _ y 
2u y ' 



( [T0| ) can be written as. 



u 



-(AV - AXu + l)z 2 + -K(2 + \u)e^ u + -KA 
3 4 2 



(A + l) 



(41) 
(42) 



wich can be studied numerically. 

Using fl2"5|), we can determine R(t) 



A ,. 



— = a0exp-0 . 



~ y 
y ~ y. 



If A ~ 0, equation @ leads to fl3|) with J 
that A; = 0. If 7^ 0, the case A ~ can be studied from: 

y Ay K . . 

y 3y 2 

With the above notation this gives 

^~ 3^ 2 + y(A+l). 
The case A = — 1 is obvious. When A > — 1, we get: 



(43) 

\K(1 + A) if we assume 



(44) 



z(t) 



(^-( A+1 )J tan 



1 (ZK 



1/2 



(A+l) (t-t ) 



or 



v(*) 



cos 



1 (ZK 



1/2 

(A + 1) ) (t - £ c ) 



+ c 



where c is a constant. We have thus: 



0(t) = exp<j3(^(A + l) 



-1/2 sin 



§(f (A + l)) 1/2 (t-to 



x 



. 7T 1 /3i^T . 

tan ' 4 + 6 llP °' 



2 cos 2 

1/2' 



|(f(A + l)) 1/2 (t-t ) 

(A + D)- 1/2 

exp(ct). (45) 



3 I ZK 
2 I 2 
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When A < — 1 we have: 



exp 



+ - tan 
2 



-3('^(A + 1) 



-1/2 sinh 



^(f(A + l)) 1/2 (t-t D 



2 cosh 2 



1/2 



sinh^^(-l-A)) (/-/„) 



^(f(A + l)) 1/2 (t-to) 

exp(ct). (46) 



4.3 A F-L universe with matter: p = up 

Let us multiply (|27| ) by and let us add (|27| ) and fl2"8|). We get: 

_) pu, + 1) + _ + _(3u, + 1) 
+ ye— (1 + - Acu0 2 ) + —K(uj + 1) = 0. 



Using (|29|) , we find: 



^-i(l + A0 2 ) + £(A 2 4 -A0 2 + 4). 



# 3(f) 



(47) 



(4* 



Then, by (^) and using the preceding notation, we can write the the 
equation: 

y 2 /(3 + uj) + 2Xujy 2 + X 2 (l+uj)y 4 \ 3yKA(uj + l) 
V J\ 2(1 + Ay 2 ) 7 + 4 1 + Xy 2 



3,, a„2 (1 + w)y — Xcuy 3 3 a„ 2 

2^ - — 1 — pi y 

4"~ 1 + Xy 2 2 - » " l + Ay 2 



+ ge^a^V^A 173 ^^.^ (49) 



which is nothing but ( j39|) when u = and p = 0; p = implies a = (see 
(H3|)). Using the conservation law V„T^ t = and equation(E9j), we find: 



= _ 3 « = (1±W (50) 

*>(! + «) R v 

An example of numerical solutions is given in Fig. 3, 4 and 5. 

In Fig. 3, we can see the typical behaviour of a Friedmann-Lemaitre uni- 
verse whose expansion is driven by exponentially harmonic maps and matter 
or radiation. This evolution know three different ages: first, the universe 
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behaves in the usual way under the presence of matter or radiation, then 
its characteristic deceleration stops and finally, when the exponentially har- 
monic map <p becomes too small, the expansion accelerates (see Eg. (|4"5D ) . But 
this is also the typical bahaviour of a universe that is dominated eventually 
by a positive cosmo logical constant or something that acts like it, such as 
a potential related to a scalar field (quintessence). Nevertheless, the careful 
reader has already noticed that this particular behaviour is only due to the 
exponentially harmonic feature of the scalar field <fi and its particular ini- 
tial conditions 4>q and 0o- This proves again that these interesting features, 
suggested among others by the observations of type la supernovae, arise nat- 
urally in the context of exponentially harmonic maps without requiring an 
adhoc mechanism of inflation. 

In Fig. 4, we have represented the evolution of the energy densities of 
matter, radiation and the scalar field coupled to it. In the very early ages of 
the universe^, the dominant contribution to the total energy is due to (f>. But 
the universe expands just as if there was only usual matter and radiation. As 
this expansion goes on, the energy density of the scalar field goes below 
the density of matter and radiation, making the expansion accelerating when 
the field <fi become tiny (see (^3[)). This would correspond to the well known 
"cosmic coincidence". Note that we have chosen the arbitrary constants 
arising in the integration of Eq.fl5U|) in order to get this peculiar model of 
expansion. The present study is purely qualitative, more work is under way 
to construct a quantitative model that could be compared, for example, with 
those of quintessence. 

In Fig. 5, we have represented the evolution of the scalar field <fi in order 
to satisfy Eq.(f4~9|) both in the presence of matter (lu = 0) and radiation 
(w = l/3). 

5 Conlusions 

The main interest of the use of exponentially harmonic maps is the fact that 
the Lagrangian ( p6|) is a generalisation of bosonic string Lagrangian (with 
only a dilatonic field) written in Einstein frame ||. Indeed, if A ~ 0, the 
Lagrangian (pE|) tends to this bosonic string Lagrangian, if A = —1 — A D , 
where A G is the small cosmo logical constant. But mathematically, this limit 
A — > is highly problematic due to the fact that unlike usual harmonic 
maps, exponentially harmonic maps have no good invariance property under 
homothetic changes of the field, as it was shown in the introduction. 

1 It is interesting to note that here there is a huge difference in the time scales of Figs. 
3 and 4, about 7 orders of magnitude. 
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An interesting open question is whether solutions of equations deriving 
from fl26|) admit duality symmetries. We are in the process of studying this 
question. 
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Figure 1: Behaviour of the expansion factor R(t) and the exponentially harmonic map 
4>(t) on an euclidean Friedmann-Lemaitre spacetime (X w ; A = and (j) > 0). 




Figure 2: Behaviour of the expansion factor R(t) and the exponentially harmonic map 
<f>(i) on an euclidean Friedmann-Lemaitre universe in the Einstein frame. (X w ; A = 0; 
e > \c\ +b anda = - {^f' 2 ). 
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1e+08 2e+08 3e+08 4e+08 5e+08 6e+08 

t 

Figure 3: Expansion factor computed from Eqs. 49 and 43 illustrating the behaviour of 
an hypothetic universe filled with matter (to = 0) and an enxponentially harmonic map. 
(<h = 1, 4>o = -5, w = 0). 
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Figure 5: Exponentially harmonic maps <f) — ±\A* solution of Eq.49, for lo = (dust- 
dominated universe) and u = 1/3 (radiation-dominated universe), ((po — 1, O = ~5 ; 
A = 0.1J. 
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